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ABSTRACT 


These  technical  memoranda  (ARL-TM-73-11  and  ARL-TM-7J-12) 
describe  two  analyses  of  ACODAC  data  tapes.  The  first 
memorandum  contains  the  initial  or  preliminary  analysis  of 
an  ACODAC  performed  at  ARL.  Data  generated  by  four  types  of 
sources --ambient  noise,  ew  source,  vibroseis  source,  and 
shot  (explosive)  source — -were  analyzed.  The  second  memorandum 
describes  the  analysis  performed  and  results  obtained  at 
ARL  in  support  of  the  Test  Director  for  the  Blake  Test  to 
evaluate  ACODAC  system  self-noise  effects  on  system  per- 
formance . 


Preceding  page  blank 


TABLE  OF  CONTENTS 


ABSTRACT 

LIST  OF  ILLUSTRATIONS 
LIST  OF  TABLES 

I.  ARL  PRELIMINARY  DATA  ANALYSIS  FROM  ACODAC  SYSTEM 
(ARL-TM-75-1L) 


A. 

B. 

C. 

D. 

E. 

F. 


G. 

H. 

I. 

J. 


Introduction 

Playback  System  and  Procedures 

Ambient  Noise 

Calibration 

Vibroseis 

cv  Sources 

1.  Spectral  Analysis 

2.  Adaptive  Spectra  Estimation 
Validation  of  Data 

Digitization  and  Analysis  of  Shot  Lata 
Intensity  Spectrum  of  Shot  Data 
Estimation  of  Signal-to-Noise  Ratio 


Page 

iii 

vii 

x 

1 

1 

1 

4 

8 

8 

12 

12 

19 

19 

21 

30 

30 


II.  ANALYSIS  OF  THE  BLAKE  TEST  ACODAC  DATA  ( AKL-TM-73-12) 

A.  Blake  Test  Analysis 

1.  Introduction 

2.  Data  Reduction  Procedures 
J.  Analysis  Constraints 

4.  Dbta  Tape  Quality 

B,  ACODAC  Calibration  Analysis 

1 . 5\irpose 

2.  Data  Seduction  Procedures 

3.  Spectral  Analysis 


1 

1 

1 

1 


3 


3 


S0. 


Renting  page  blank 


v 


'•'■if- -• 


TABLE  OP  CONTENTS  (Coat’d) 


Page 

C.  Ambient  Noise  Analysis 

13 

1. 

Data  Analyzed 

13 

2. 

Data  Reduction  Procedures 

13 

3. 

Spectra 

13 

D.  Strumming  Analysis 

13 

1. 

Purpose 

15 

2. 

Data  Reduction  Procedures 

15 

3. 

Results 

26 

E.  Shot 

. Data  Analysis 

26 

1. 

Introduction 

26 

2. 

Multipath  Structure  of  Signals 

32 

3. 

Signal -to-Noise  Ratio 

32 

4. 

Spectrum 

48 

APPENDIX 

A 

31 

APPENDIX 

R 

61 

LIS?  OF  ILLUSTRATIONS 


!•  ARL  PRELIMINARY  DATA  ANALYSIS  FROiM  ACODAC  SYSTEM 

( ARL-TM-73-H) 


Figure 


Title 


Page 


1 

3 

4 

3 

6 

7 

8 

9 

10 

11 

12 

15 
1* 
Is? 

16 
17 
15 
19 

3Q 


<zd' 


S3 


Applied  Research  Laboratories  CDC  3200  Computer 
Configuration 

ARL  Data  Reduction  Facility 

Ambient  Noise  Spectra 

Ambient  Noise  Spectra 

Ambient  Noise  Spectra 

Ambient  Noise  Time  Varying  Spectra 

ACODAC  Calibration  Signal  Spectra 

Vibroseis  Source  Time  Varying  Spectra 

Vibroseis  Source  Time  Varyiiig  Spectra 

e«  Source  Spectra 

cw  Source  Spectra 

cv  Source  Spectra 

cv  Source  Spectra 

cw  Source  Time  Varying  Spectra 

cv  Source  Adaptive  Spectral  t tt  "ration 

Test  for  HottOiceneity  Ambient  »k-*  • cv  Source 

Representative  Ceopu&er  Piet#  of  Digital  Shot  Data 

Representative  Computer  Hot#  of  Digital  Shot  Data 

Comparison  of  Structure?  ef  hscelvei  Shot  Signal# 
frees  beginning  and  End  of  Data  Set# 

Detailed  ?<*spwi#on  of  Two  Primry  Arrival#  fre®:- 
Shot  Sources 

Effect  of  deduction  ir,  Saapl*  Rate  on  Preservation 
of  ifaveTcrst  Structure 

C«e&arlso«  of  the  Vaveforis  Structures  of  the  Second 
Arrivals,  of  the  First  Four  Shots 

Intensity  Spectrum  of  the  First  Arrival  cf  the  F**urth 
Sh©t  Sequence 


2 

3 

5 

6 
7 
9 

10 

11 

13 

14 

15 

16 

17 

18 

?o 
( ^ 
23 


31 


LIST  OF  ILLUSTRATIONS  (Cont'd) 


Figure 

Title 

Page 

2b 

Estimation  of  Shot  Signal-to-Noise  Ratio,  S/N 

55 

25 

Signal-to -Noise  Ratio  Estimates  as  Functions  of 
Averaging  Time  of  Noise 

54 

II.  ANALYSIS  OF  THE  BLAKE  TEST  ACODAC  DATA 
(ARL-TM-73-12) 

Figure 

Title 

Page 

1 

Calibration 

Spectra  for  WHOI  Array  - Hydrophone  l 

4 

2 

Calibration 

Spectra  for  WHOI  Array  - Hydrophone  2 

5 

5 

Calibration 

Spectra  for  WHOI  Array  - Hydrophone  5 

6 

h 

Calibration 

Spectra  for  WHOI  Array  - Hydrophone  4 

7 

5 

Calibration 

Spectra  for  WHOI  Array  - Hydrophone  5 

8 

6 

Calibration 

Spectra  for  WHOI  A ray  - Hydrophone  6 

9 

7 

Calibration 

Spectra  for  UM  Array  - Hydrophone  6 

10 

8 

Ambient  Noise  Spectra  for  WHOI  Array  - Hydrophone  1 

16 

9 

1/5  Octave  Spectra  of  Ambient  Noise  for  WHOI  Array  - 
Hydrophone  1 

17 

10 

Ambient  Noise  Spectra  for  WHOI  Array  - Hydrophone  2 

18 

11 

Ambient  Noise  Spectra  for  WHOI  Array  - Hydrophone  5 

19 

12 

Ambient  Noise  Spectra  for  WHOI  Array  - Hydrophone  4 

20 

15 

Ambient  Noise  Spectra  for  WHOI  Array  - Hydrophone  5 

21 

14 

'mbient  Noise  Spectra  for  WHOI  Array  - Hydrophone  6 

22 

15 

Ambient  Noise  Spectra  for  UM  Array  - Hydrophone  6 

25 

16 

Time  Varying 

Spectra  for  WHOI  Hydrophone  1 

26 

17 

Time  Varying 

Spectra  for  WHOI  Hydrophone  2 

27 

18 

Time  Varying 

Spectra  for  WHOI  Hydrophone  5 

28 

19 

Time  Varying 

Spectra  for  WHOI  Hydrophone  4 

29 

20 

Time  Varying 

Spectra  for  WHOI  Hydrophone  5 

30 

viii 


LIST  OF  ILLUSTRATIONS  (Cont'd) 


Title 


2i 

Time  Varying  Spectra  for  WHOI  Hydrophone  6 

31 

22 

Amplitude  versus  Time,  Digitized  Record  4-Shot 
Sequence,  10  mile  Range,  300  ft  Shot  Depth, 
WHOI  Hydrophone  1,  First  Half 

33 

23 

Amplitude  versus  Time,  Digitized  Record  4-Shot 
Sequence,  10  mile  Range,  300  ft  Shot  Depth, 
WHOI  Hydrophone  1,  First  Half 

34 

24 

Amplitude  versus  Time  Record  of  Two  Shots  Showing 
Multipath  Arrivals,  10  mile  Range,  WHOI  Hydrophone  1 

35 

25 

Estimation  of  Shot  Signal-to-Noise  Ration  (S/N) 

36 

26 

Shot  Signal  and  Signal-to-Noise  Ratio 

39 

27 

Signal-to-Noise  Ratio  of  Shot  Data  as  a Function  of 
WHOI  Hydrophone  Number 

4o 

28 

Signal-to-Noise  Ratio  of  Shot  Data  as  a Function  of 
WHOI  Hydrophone  Number 

4i 

29 

Signal-to-Noise  Ratio  of  Shot  Data  as  a Function  of 
WHOI  Hydrophone  Number 

42 

30 

Signal-to-Noise  Ratio  of  Shot  Data  as  a Function  of 
WHOI  Hydrophone  Number 

43 

31 

Variation  of  S/N  with  Range  of  Shot 

44 

32 

Variation  of  S/N  with  Range  of  Shot 

45 

33 

Variation  of  S/N  with  R'inge  of  Shot 

46 

3*+ 

Comparison  of  WHOI  and  Miami  Hydrophone  6 Shot  S/K 

4? 

35 

Shot  Spectrum 

4y 

36 

Shot  Spectrum 

50 

ix 


LIST  OP  TABLES 

II.  AKAL1SIS  OP  THE  BLAKE  TEST  ACODAC  BATA 
(ABL-W-73-12) 


Title 

Signal  Level  Check  for  the  ACODAC  Tapes 
Calibration  Spectra  Labels 
Calibration  Spectral  Levels  In  decibels 
Calibration  Spectra  Labels 

Aobient  Hoise,  Relative  Spectral  Levels  in  decibels 
Shot  Pattern 

Signal  -to-Soise  Ratios : Shot  Bata 


I.  ML  PRELIMINARY  DATA  ANALYSIS  FROM  ACODAC  SYSTEM 

(ARL-TM-75-11) 


A.  Introduction 

The  duplicate  data  tape  from  the  ACODAC  System  used  during  the 
Church  Gabbro  Exercise  was  received  at  Applied  Research  Laboratories 
(AEL),  The  University  of  Texas  at  Austin,  early  on  11  June  1973.  A 
preliminary  analysis  was  conducted  to  demonstrate  ARL’s  capabilities  to 
handle  data  analysis  of  this  type.  The  analysis  was  constrained  by 
limited  time,  inability  to  read  time  code  information  and  amplifier  gain 
changes,  and  no  log  or. record  of  the  content  of  the  data  tape. 

After  examining  the  tape,  the  decision  was  made  to  isolate  and 
analyte  four  different  data  segments  containing  ambient  noise,  continuous 
wave  (cw)  sources,  shots,  and  the  vibroseis  source.  Time  permitted  only 
the  analysis  of  data  from  hydrophone  number  3. 

S.  Playback  System  and  Procedures 

Figures  1 and  2 are  block  diagrams  of  the  AEL  computer  and  data 
.-eduction  facility. 

Since  do  log  was  provided  with  chc  data  tape,  it  was  necessary  to 
search  the  tape  to  locate  four  different  types  of  data  for  analysis. 

The  analog  data,  was  played  bad  at  a speedup  cf  80:1  to  locate  different 
data  serpents  for  analysis.  Simultaneous  visual  observations  of  the 
output  of  hydrophones  2 and  3,  with  an  audio  presentation  of  the  out- 
put of  hydrophone  3,  enable a the  operator  to  locate  the  data  segments. 


A Spectral  Dynamics  analyzer  (500  line)  was  then  used  to  examine  the 
spectral  characteristics  of  these  data  segments  to  identify  ambient 
noise,  cw  sources,  shots,  and  vibroseis  source. 

These  data  segments  were  passed  through  a 10  to  300  Hz  (real  time) 
filter  to  prevent  aliasing  and  to  eliminate  extraneous  low  frequency 
playback  system  noise  outside  the  frequency  band  of  interest.  The  data 
were  digitized  at  a 600  Hz  (real  time)  rate  in  8192  sample  blocks  and 
stored  on  digital  tape.  The  600  Hz  sync  signal  for  the  A/D  converter 
was  obtained  from  the  12th  harmonic  of  the  50  Hz  carrier  on  the  time 
code  channel,  thus  eliminating  wow  and  flutter  induced  by  the  recording/ 
playback  processes. 

Additional  data  were  processed  via  the  Spectral  Dynamics  analyzer. 
The  spectral  lines  computed  by  the  analyzer  were  output  to  the  digital 
computer  and  stored  on  digital  tape. 

The  location  in  time,  relative  to  the  start  of  the  data  tane , of 
the  data  segments  analyzed  are  approximate  since  no  tape  log  was  avail- 
able. The  estimated  times  were  taken  f^om  the  reel  turu  count  of  the 
playback  recorder,  which  is  not  a linear  system. 


C 


Ambient  Noise 


Figures  3?  4,  and.  5 are  relative  power  spectra  of  three  consecutive 
segments  of  ambient  noise  data  with  no  amplifier  gain  changes  observed 
visually.  These  data  segments  occurred  toward  the  end  of  the  first 
recording  day.  The  data  were  digitized  using  the  phase  locked  600  Hz 
signal,  spectra  were  computed  with  in  FIT  using  a Hanning  window,  and 
the  spectra  were  averaged  to  produce  stable  spectral  estimates.  Each 
spectrum  covers  a data  block  of  13*7  pec  (8192  samples).  The  frequency 
resolution  is  0.07  Hz.  Each  graph  represents  an  average  of  consecu- 
tive soectra  or  a real  time  interval  of  11.4  min.  The  general  spectra 
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shape  remains  consistent  throughout  the  analysis  sequence.  The  spectral 
lines  between  25  and  75  Hz  indicate  the  presence  of  a ship  in  the  area. 
Figure  5,  the  beginning  of  the  analysis  sequence,  shows  a spectral 
increase  of  approximately  4 dB  in  the  frequency  range  of  100  to  125  Hz. 

Figure  6 is  a time  varying  power  spectra  of  a different  data 
segment  of  ambient  noise  occurring  near  the  end  of  the  first  recording 
day.  The  relative  power  spectra  were  computed  by  the  Spectral  Dynamics 
analyzer  with  a resolution  of  O.75  Hz.  The  display  shows  the  frequency 
markers  from  the  analyzer  for  reference,  the  frequency  spectra  computed 
for  the  calibration  signal  recorded  on  the  iCODAC  tape,  and  the  spectra 
of  the  ambient  noise.  Each  time  trace  represents  an  average  of  four 
power  spectra  computed  with  a 2 sec  data  block.  The  lack  of  low  frequency 
lines  indicates  no  shipping  present. 

D.  Calibration 


Figure  7 is  a spectral  analysis  of  one  of  the  calibration  signals 
occurring  every  6 h on  the  ACQBAC  recording  system.  This  data  segment 
occurred  approximately  1 1/5  days  into  the  recording.  The  data  were 
digitized  at  a 600  Hz  rote  in  8192  samples /sequence  for  1 8 sequences  or 
in  data  blocks  of  1J.7  sec.  The  frequency  resolution  in  0.07  Hz . The 
total  length  of  the  data  sequences  is  4.2  oin.  A spectns*;  is  computed 
for  each  sequence,  and  the  sequences  (18)  are  averaged  t«  achieve  stable 
estimates.  The  calibration  signal  appeared  to  be  performing  its  proper 
function. 


E.  Vibroseit 


Figure  8 is  a display  of  the  tie*  varying  relative  linc&r  power 
spectra  of  a sequence  of  vi broads  data  occurring  approximately  9 days 
from  start  of  recording.  The  relative  power  spectre  were  raesputed  with 
a Spectral  Dynamics  real  tisse  analyzer  that  produces  spectral  lines 
with  0.75  liz  resolution.  Each  spectrum  is  cocputed  for  a ? sec  data 
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block.  The  analyzer  is  set  to  average  four  consecutive  spectra, 
representing  the  analysis  of  8 sec  of  data,  and  output  to  the  digital 
computer.  These  spectra  were  again  averaged  to  represent  power  spectra 
for  to  sec  data  blocks.  The  display  is  linear  power  versus  frequency 
versus  relative  time . Frequency  markers  from  the  analyzer  are  displayed 
for  reference.  Figure  9 is  a similar  display,  but  with  log  power. 

F.  cw  Sources 

A data  segment,  containing  an  apparent  cw  source  of  9°  Hz,  was 
analyzed.  The  data  segment  occurred  2 1/2  days  into  the  recording  period. 
Two  different  types  of  analysis  were  performed:  high  resolution  spectral 

analysis  and  adaptive  spectral  estimation  to  determine  the  spectral 
content  of  the  data  segment. 


High  resolution  power  spectra  of  a cw  source  in  the  ambient 
noise  background  is  shown  in  Fig.  10,  The  spectra  are  computed  with  an 
FFT  algorithm  and  averaged  in  blocks  of  50  (11.4  min).  The  frequency 
resolution  is  0.07  Kz. 


Spectral  lines  between  10  and  75  Hz  indicate  the  presence  of 
shipping.  The  line  at  9^.7  H2  is  assumed  to  be  the  cw  source.  The 
side  lobe  structure  of  the  spectral  line  indicates  an  amplitude  modulation. 


Three  consecutive  sets  of  data  were  analyzed  in  the  same 
manner  to  observe  if  the  lines  were  present  over  a longer  time  interval. 
Figures  11,  12,  and  13  show  how  the  lines  vary  over  a 33  min  data  block. 

Time  varying  power  spectra  of  the  ambient  noise  plus  cw  signal 
are  shown  in  Fig.  14.  T>ie  spectra  are  computed  via  an  FFT  using  a 
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Gaussian  window  with  30  equivalent  degrees  of  freedom.  No  ensemble 
averaging  is  used.  The  basic  frequencj  resolution  is  G.29  Hz.  Each 
line  is  an  FFT  on  the  first  2048  points  from  every  3rd  sequence  on 
digital  tape  1604.  The  cw  tone  at  90*7  Hz  is  evident  and  it  is  changing 
in  emplitude  with  time. 

2.  Adaptive  Spectra  Estimation 

Figure  10  shows  definite  spectral  lines.  In  an  effort  to 
check  on  the  stability  of  these  lines,  a modification  of  the  maximum 
entropy  method  (MEM)  of  computing  spectra  was  used  on  the  first  4l 
records  of  digital  tape  1604.  The  basic  program  forms  an  autocovariance 
matrix,  inverts  it,  and  solves  for  the  roots  of  a polynomial,  which 
indicate  resonance  frequencies  in  the  data.  The  primary  advantage  of 
the  method  is  that  it  allows  the  use  of  a small  amount  of  data  to  obtain 
a high  resolution  spectrum  or  the  resonance  frequencies  in  the  case  of 
this  modified  version.  The  amplitude  of  spectral  lines  is  ignored  by 
the  program,  but  they  may  be  found  later. 

Figure  15  shows  the  results  of  a time  versus  frequency  display. 
The  program  is  cocrraanded  to  search  for  l4  roots;  and  up  to  l4  asterisks, 
which  represent  the  frequencies  found  at  that  time,  are  displayed  on  each 
line.  The  90  Hz  line  is  strong  and  consistent;  the  other  lines  appear 
to  fade  in  and  out.  The  method  is  dependent  upon  signal-to-noise  ratio 
(S/N),  and  as  the  signal  fades  the  lines  will  appear  to  shift  frequency. 
The  conclusion  is  that  the  lines  are  present  on  the  average,  and  they 
fade  in  and  out  in  magnitude. 

G.  Validation  of  Data 


A useful  measure  of  the  changing  statistics  of  noise  data  is  to 
test  for  homogeneity.  That  is,  if  two  data  sets  have  the  same  proba- 
bility density,  then  they  are  homogeneous.  In  this  case,  a long  digital 
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tape  subdivided  into  records  of  8192  samples  each  was  examined.  This 
data  segment  occurred  21/2  days  into. the  tape  and  contains  ambient  noise 
plus  a cw  signal.  The  10th  record  was  chosen  as  a reference  and  was 
compared  to  each  of  the  other  138  records . One  hundred  samples  were 
obtained  from  each  record  by  picking  every  10th  sample  from  the  first 
1000  samples.  The  Kolmogorov-Smirnov  (K-S)  statistic  was  computed  for 
each  comparison  and  the  result  is  shown  in  Fig.  16.  The  conclusion  is 
that  after  a 2 h span  of  time  the  statistics  change  sufficiently  for 
the  K-S  test  to  consistently  reject  the  hypothesis  that  the  data  are 
homogeneous.  After  4 h the  data  are  again  homogeneous  with  respect 
to  the  10th  record.  This  test  should  be  repeated  using  different  records 
as  references  for  a complete  analysis. 

This  brief  analysis  shows  that  for  this  environmental  situation  it 
would  not  be  appropriate  to  average  for  more  than  a 2 h span  of  time 
to  obtain  a statistically  reliable  estimate  of  power  spectra. 

H.  Digitization  and  Analysis  of  Shot  Data 


To  obtain  a suitable  data  base  that  would  facilitate  a preliminary 
investigation  of  the  shot  signals,  a total  of  36  waveforms  from  explosive 
sources  were  digitized.  These  data  were  approximately  5 l/2  days  into  the 
tape  and  were  generated  with  a "2- shot"  explosion  seqv.ence.  Each 
received  signal,  from  hydrophone  J,  was  gated  in  time  by  choosing  an 
appropriate  time  delay  following  the  first  shot  in  each  pair.  The  sample 
rate  was  the  phase  locked  600  Hz.  Figure  17  shows  computer  plots  of  the 
first  five  echoes  of  the  56-echo  sequence.  The  location  of  the  shot 
signals  are  evident  although  only  every  50th  sample  is  plotted  here.  The 
reason  for  choosing  a to  sec  window  to  select  the  data  was  to  overcome 
instability  of  the  triggering  method.  More  detailed  plots  of  the  data 
are  shown  ir.  Fig.  18.  Here  every  5th  sample  has  been  plotted. 


Kolmogorov  - Smirnov 
Statistic 
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A general  structure  of  the  shot  data  is  shown  in  Pig.  19  where  three 
consecutive  signals  free  the  beginning  of  the  data  set  (of  36  shots)  and 
two  consecutive  signals  from  the  end  of  the  data  set  are  plotted.  The 
first  three  signals  are  similar  in  structure  as  are  the  last  two  signals. 
However,  a structural  difference,  with  respect  to  the  number  and  arrival 
times  of  multipath  signals,  does  exist  between  the  set  of  the  three 
earlier  shots  and  the  set  of  two  later  shots.  This  is  to  be  expected 
since  the  geometry  of  the  shot  experiment  changed  over  the  duration  of 
the  measurements,  and  accordingly,  the  multipath  structure  would  also  change. 
Within  a few  consecutive  shots,  however,  the  multipath  structure  appears 
to  be  relatively  stable. 

Figure  20  provides  a comparison  for  the  waveforms  of  the  primary 
arrivals  of  two  (Ho.  1 and  5)  shot  signals.  The  two  waveforms  are  very 
similar  ir  shape , particularly  with  respect  to  the  location  and  amplitudes 
of  some  of  the  larger  peaks . Some  differences  between  these  two  waveforms 
can  be  discerned,  and  these  could  be  attributed  to  either  differences  in 
the  generation  of  the  shots  or  changes  in  the  propagation  environment. 

To  obtain  a preliminary  estimate  of  the  required  sample  rate  on  the 
shot  data,  a signal  was  plotted  with  various  sample  rates.  It  was  con- 
cluded that  at  least  600  He  is  necessary  and  very  likely  sufficient  for 
the  preservation  of  the  information  in  the  signal  waveforms.  Figure  21 
illustrates  the  effect  of  reducing  the  sample  rate  on  a typical  signal. 
Misleading  results  will  be  obtained  with  a $00  fte  sample  rate  since  much 
of  the  waveform  structure  is  lost  when  the  sample  rate  is  changed  from 
600  Hz  to  303  Hz. 

Using  Fig.  22  one  can  compare  the  shot -to-shot  consistency  of  the 
second  arrivals.  The  consistency  here  does  not  appear  to  be  as  high  as 
coat  of  the  primary  arrival  (see  Fig.  20).  Also,  the  second  arrivals 
are  more  spread  in  time,  which  indicates  that  this  particular  multipath 
component  introduces  a time  smear  into  the  acoustic  propagation. 
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FIGURE  20 

DETAILED  COMPARISON  OF  TWO  PRIMARY 
ARRIVALS  FROM  SHOT  SOURCES 
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FIGURE  21 

EFFECT  OF  REDUCTION  IK  SAMPLE  PATE  ON 
PRESERVATION  OF  WAVEFORM  STRUCTURE 
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I.  Intensity  Spectnm  of  Shot  Data 


A very  limited  attempt  was  made  to  estimate  the  intensity  spectrum 
of  the  received  shot  data.  A sample  spectrum  (2,34  Hz  resolution)  of 
the  first  arrival  of  the  fourth  shot  is  shown  in  Fig.  23.  Generally, 
the  spectrum  decays  from  its  maximum  near  zero  frequency  to  a minimum 
at  300  Hz,  the  extent  of  the  receiver  filters.  This  estimate  has  many 
fluctuations  present,  some  of  which  are  undoubtedly  due  to  statistical 
inaccuracies  in  the  estimate.  Time  permitting,  the  estimate  would  be 
improved  by  several  techniques  such  as  averaging  several  spectra  or  using 
appropriate  smoothing  functions. 

Before  proceeding  with  the  computation  of  the  spectra,  however,  the 
shot  data  should  be  more  accurately  described  with  respect  to  its 
deterministic  or  random  structure.  If  the  shot  data  were  entirely  determi- 
nistic and  repeated  consistently  over  the  duration  of  the  experiment,  then 
al’  of  the  spectral  information  could,  in  principle,  be  obtained  from  one 
sample  function.  However,  the  shot  data  do  appear  to  have  random  com- 
ponents xnd  particularly  so  in  the  2nd,  3rd,  4th,  etc.,  arrivals.  When 
this  is  the  case  then  appropriate  ensemble  techniques  must  be  employed 
to  obtain  a statistically  reasonable  estimate  of  the  intensity  spectra  of 
the  data.  This  technique  should  be  used  carefully,  however,  since  the 
"time"  spectrum  of  the  shot  data  may  change  during  the  course  of  an 
experiment  due  to  a change  in  the  experimental  geometry  and  propagation 
conditions.  In  this  case,  ensemble  spectra  could  be  misleading  since 
it  would  be  derived  from  data  the,  were  essentially  inhomogeneous.  This 
problem  can  only  be  solved  after  a preliminary  investigation  of  the  sta- 
tistical properties  of  the  received  shot  data. 

J.  Estimation  of  Signai-to-Kolss  Ratio 

One  of  the  more  important  quantities  id  the  S/S5  of  any  given  set 
shot  data.  In  this  analysis  the  S/E  of  a typica:  signal  was  estimated 
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FIGURE  23 

INTENSITY  SPECTRUM  OF  THE  FIRST  ARRIVAL  OF 
THE  FOURTH  SHOT  SEQUENCE 
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According  to  two  definition*  of  S/N  and  several  averaging  times  used  in 
the  estimation  of  the  noise  background.  The  problem  is  illustrated  in 
Pig.  24,  In  thir  case  the  signal  is  relatively  easy  to  locate,  a situa- 
tion that  likely  will  not  always  occur.  The  epoch  of  the  signal  is 
noted  as  well  as  a particular  window  used  for  the  estimate  of  the  noise 
power  or  intensity.  Ihe  two  S/N's  estimated  are  computed  using  the  peak 
of  the  signal  and  the  average  energy  of  the  signal  (see  Pig.  25).  The 
stability  of  the  estimate  should  improve  as  the  averaging  time  of  the 
noise  estimate  improves.  This  is  illustrated  very  clearly  in  Fig.  25 
where  It  can  be  seen  that  the  S/N  estimate  approaches  a stable  value  when 
an  averaging  time  of  0.5  sec  or  greater  is  used.  To  determine  useful 
S/N  values,  analysis  of  the  effect  of  averaging  time  should  be  performed  to 
obtain  higher  confidence  in  the  results  as  well  as  a savings  in  processing 
time. 
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FIGURE  25 

SIGNAL  TO  NOISE  RATIO  ESTIMATES  AS 
FUNCTIONS  OF  AVERAGING  TIME  OF  NOISE 
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II.  ANALYSIS  OF  THE  BLAKE  TEST  ACODAC  DATA 
(AEL-TM-73-12) 


A.  Blake  Test  Analysis 

1.  Introduction 


The  Blake  Test  was  conducted  at  sea  in  a deep  water  environment 
for  the  purpose  of  evaluating  the  ACODAC  system  self-noise  effects  on 
system  performance.  Data  tapes  from  the  Woods  Hole  Oceanographic  Institute 
(WHOl)  and  the  University  of  Miami  (UM)  ACODAC  systems  were  duplicated 
and  furnished  to  AKL  for  analysis  in  support  of  the  Test  Director  (UM). 

The  UM  array  parted  early  into  the  test  leaving  only  hydrophone  channel  6 
recording  daca. 

2.  Data  Reduction  Procedures 


Data  reduction  was  accomplished  using  the  hardware/software 
systems  currently  in  operation  at  ARL.  These  systems  are  described  in 
Section  I (ARL-TM- 73-11.) . The  variances  on  these  basic  procedures  are 
outlined  in  this  report  with  the  description  of  each  type  of  analysis 
performed. 

3.  Analysis  Constraints 

The  data  analysis  was  somewhat  constrained  by  the  nature  of  the 
test,  which  included  a tight  time  schedule  for  the  completion  of  the 
analysis.  The  data  tapes  were  received  at  ARL  on  20  June  1973  and  the 
results  were  presented  at  a meeting  with  the  Test  Director  at  WHOI  on 
25  June  1973- 
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Other  constraints  that  slowed  the  analysis  were  the  lack  of  a 
time  code  reader  and  test  event  log.  The  ACODAC  tapes  were  edited  and 
consistent  events  were  located  on  each  tape.  For  reference  in  coordinating 
this  analysis  with  other  analyses,  the  time  code  was  digitized  by 
quadrature  sampling  and  decoded  with  the  digital  computer.  The  ACODAC 
amplifier  gains  were  read  with  the  same  procedure. 

4.  Data  Tape  Quality 

The  duplicate  ACODAC  analog  tapes  were  examined  for  "quality" 
prior  to  the  analysis.  The  signal  levels  on  the  Blake  Test  tapes  were 
compared  to  the  Church  Gabbro  tape  previously  analyzed  (Section  i). 

Table  I gives  the  signal  level  on  the  tape,  using  identical  playback 
electronics,  for  the  Church  Gabbro,  Blake  WHOI,  and  Blake  UM  tapes.  The 
peak  calibration  signals  for  the  Blake  Test  tapes  were  10  to  20  dB  below 
comparable  signal  levels  on  the  Church  Gabbro  tape.  The  ARL  playback 
system  noise,  using  blank  tape  in  the  recorder,  is  also  given  in  Table  I. 

TABLE  I 

SIGNAL  LEVEL  CHECK  FOR 
THE  ACODAC  TAPES 
Signal  Levels  (dB  re  1 Vrms) 


Hydrophone 

Channel 

Playback* 

system 

noise 

Church 

Gabbro 

WHOI 
( Blake ) 

UM 

Blake 

1 

-1*9 

-11 

-21 

- 

2 

-45 

-9 

-21 

- 

3 

-49 

-11 

-33 

- 

1* 

-47 

-10 

-27 

- 

5 

-47 

-15 

-43 

- 

6 

-47 

-9 

i 

ro 

vo 

-43 

•Playback  system  noise  was  determined  by  using  blank  tape. 
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B.  ACODAC  Calibration  Analysis 


1.  Purpose 

The  analysis  time  period  for  the  Blake  Test  data  did  not  permit 
an  extensive  study  of  the  stability  of  the  ACODAC  system  calibration.  To 
establish  a "quality  check"  on  the  data  for  each  hydrophone  to  be  analyzed, 
the  calibration  signal  recorded  immediately  prior  to  the  ambient  noise 
data  was  examined.  This  calibration  occurred  between  the  50  and  100  mile 
shot  sequences. 

2.  Data  Reduction  Procedures 


The  time  compressed  (40:l)  analog  data  was  bandpass  filtered 
(10  to  300  Hz)  and  digitized.  The  A/D  converter  was  synced  to  the  12th 
harmonic  (600  Hz)  of  the  50  Hz  carrier  of  the  time  code  to  reduce  tape 
recorder  inaccuracies.  The  spectra  were  computed  with  a FFT  algorithm  to 
a resolution  of  O.O73  Hz.  The  data  block  for  the  FFT  is  13.65  sec.  A 
Hanning  window  was  applied  to  each  spectrum  for  smoothing.  The  spectra 
shown  in  Figs.  1 through  7 represent  an  ensemble  average  of  2.3  min.  The 
spectra  have  been  corrected  for  the  playback  amplifier  gain  but  not  for 
the  ACODAC  amplifier  gain.  Table  II  is  an  explanation  of  the  labels 
for  the  spectra. 

3.  Spectral  Analysis 

Table  III  is  a summary  of  the  peak  and  noise  spectral  levels 
of  the  calibration  signal  for  the  WH0I  and  UM  arrays.  For  the  WH0I  array, 
hydrophone  channels  1,  2,  and  4 show  a consistent  peak  signal-to-noise 
ratio  (S/N)  of  at  least  50  dB  at  50  Hz.  Hydrophone  channels  3,  5.  and 
6 of  the  WH0I  array  indicate  a lower  dynamic  range.  For  all  hydrophone 
channels  the  200  Hz  calibration  signal  is  approximately  10  dB  lower  than 
the  50  Hz  calibration  signal.  All  of  the  channels  contain  harmonics  of 
the  calibration  signals  of  significant  levels. 
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TABLE  II 


CALIBRATION  SPECTRA  LABELS 


LABEL 

SYSTEM 


BXPLA3ATI0 


5 WHOI  Array 

4 


DAY 

HOUR 

MIN 


Time  Code 


CHANNEL  Hydrophone  number 

(top  hydrophone  = l) 

AMPGAIN  ACODAC  Ansplifier  Gain  as  Decoded 


SCALE  Playback  Amplifier  Gain 


SAHPH 

FREQR 

S&8QTH 


BLKTIMS 


Saspling  rate*  He 
Frequency  resolution*  Hs 
Smoothing  window 
-i  o Hanning  window 
.3311  a III  octave  filter 

Data  block  length*  sc-c 


AVGTIHE  Ensessble  averaging  ties,  fein- 


TABLE  III 


CALXBRATIOH  SPECTRAL  LEVELS  IB  DECIBELS 
WHOI  Array 


Hydrophone 


Butter 

Boise  (10-3QG) 

22. 

1 

-38 

♦17 

2 

-38 

♦IT 

3 

-42 

*5 

4 

-39 

♦11 

5 

—42 

-7 

6 

-39 

♦1 

100 

125 

W<|| 

150 

205 

?io 

-21 

-18 

-21 

♦5 

-22 

-36 

-18 

-18 

♦1 

-34 

-15 

-28 

-21 

_■? 

-30 

-13 

-23 

-22 

♦3 

-32 

-24 

-39 

-27 

-16 

-31 

-9 

-32 

-12 

-9 

-20 

UK  Array 


n 
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Hydrophone  channel  6 of  the  UM  array  showed  a peak  P/N  of 
33  dB  at  50  Hz.  The  200  Hz  cal  signal  is  13  dB  below  the  level  of  the 
30  Hz  cal  signal. 

The  broadening  of  the  cal  spectral  lines  at  50  and  200  Hz  may 
have  been  caused  by  feed  over  of  the  calibration  signal  into  the  time 
code  track  thereby  causing  the  A/D  sync  signal  to  vary. 

C.  Ambient  Holse  Analysis 

1.  Data  Analyzed 


The  data  chosen  for  the  ambient  noise  analysis  were  recorded  between 
the  50  and  100  mile  shot  sequence.  The  decoded  time  code  read  day  3, 
hour  0,  minutes  S to  46  for  the  WHO!  array  and  day  2,  hour  0,  minutes  9 
to  46  for  the  131  array. 

2.  Bata  Reduction  Procedures 


The  analog  data  on  the  A2QDAC  tapes  were  time  compressed 
(40:1)  in  playback,  bandpass  filtered  (10  to  500  He),  amplified,  and 
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digitised.  The  A/D  converter  word  is  11  bits  plus  sign  and  was  synced 
to  the  12th  harmonic  (600  Ms)  of  the  5°  B*  carrier  of  the  time  code  to 
siniaite  record/' reproduce  errors.  The  spectra  were  computed  with  a 
fft  algorithm  with  a resolution  of  O.OJJ  Us.  Each  data  block  analysed 
was  13.0  sec.  A Hanning  wi{j-f>v  was  used  for  smoothing.  The  spectra 
were  ensemble  averaged  for  11.4  « an.  Each  spectra  is  corrected  for  the 
playback  amplifier  gain  but  not  for  the  AC03AC  system  gala.  Table  tv  is 
an  'nation  of  the  labels  for  the  spectra. 
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TABLE  IV 


CALIBRATION  SPECTRA  LABELS 


LABEL 

EXPLANATION 

SYSTEM 

5 WHOI  Array 
4 UM  Array 

DAY  I 

HOUR) 

Time  Code 

MIN  J 

CHANNEL 

Hydrophone  number 
(top  hydrophone  = l) 

AMPGAIN 

ACODAC  Amplifier  Gain  as  Decoded 

SCALE 

Playba-'k  Amplifier  Gain 

SAMPR 

Sanpling  rate,  Hz 

FREQR 

Frequency  resolution,  Hz 

SMOOTH 

Smoothing  window 

-1  = Hanning  window 
•3333  3 1/3  octave  filter 


3LKTIME  Data  block  length,  sec 

AVGTIME  Ensemble  averaging  time,  min. 
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3.  Spectra 


Figures  8 through  15  give  the  ambient  noise  spectra  for  11.4  min 
of  data  as  a function  of  depth.  The  spectral  lines  at  frequencies  less 
than  50  Hz  are  assumed  to  be  due  to  ship  traffic.  Table  V gives  the 
relative  spectral  levels  for  the  peak  near  50  Hz  and  the  system  noise, 
both  record  and  playback,  at  300  Hz.  These  spectral  levels  and  the 
differences  between  the  50  and  300  Hz  levels  indicate  the  errors  resulting 
from  the  improper  operation  of  the  different  hydrophone  channels  of  the 
ACODAC  systems  and  the  analog  tape  duplication. 

D.  Strumming  Analysis 

1,  Purpose 

The  Blake  Test  was  conducted  to  evaluate  the  self -noise  or 
strunming  effects  on  the  performance  degradation  of  different  hydrophone 
array  configurations . The  WHOI  array  is  the  taut  array  type,  and  the  DM 
array  is  the  compliant  array  type.  The  UM  array  parted  very  early  in  the 
teRt;  therefore  the  only  data  presented  is  for  the  WHOI  or  taut  array 
configuration.  Spectra  of  ambient  noise  data  for  the  frequency  range  of 
0 to  6.25  Hz  were  used  to  illustrate  the  effects  of  strumming.  The  data 
segment  was  recorded  beginning  at  time  code  day  3>  hour  0,  minute  8. 

2.  Data  Reduction  Procedures 


The  analog  ACODAC  tape  was  time  compressed  (320:l)  in  playback, 
bandpass  filtered  (l  Hz  to  100  Hz),  and  amplified.  This  analog  signal 
was  input  to  a Spectral  Dynamics  SD301C  Real  Time  Analyzer  (RTA),  which 
output  500  spectral  lines  with  a frequency  resolution  of  0.19  Hz  into  the 
CDC  3200  digital  computer.  The  ensemble  averager  of  the  RTA  was  used  to 
average  the  spectra  for  64  sec.  The  reduction  procedure  consisted  of 
first  digitizing  the  frequency  markers  from  the  RTA  for  reference,  and 
then  digitizing  the  spectra  of  a data  segment  for  each  hydrophone  channel. 
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The  spectra  are  relative  in  magnitude  on  a channel-to-channel  basis  since 
the  ACODAC  system  gain  corrections  have  not  been  applied.  Playback  system 
noise  spectra  levels  are  shown  for  hydrophone  channels  5 and  6. 

3.  Results 


Figures  1 6 through  21  are  time  varying  spectra  of  the  ambient 
noise  field.  These  spectra  illustrate  the  large  levels  below  6 to 
25  Hz  due  to  the  strunming  or  self '■noise  of  the  array.  The  effects  of  the 
large  amplitude  low  frequency  signals  in  depressing  the  gains  of  the 
ACODAC  amplifiers  and  reducing  the  dynamic  range  of  the  system  are  illus- 
trated in  the  spectra.  The  effects  of  the  strunming  vary  from  hydrophone 
channel  to  hydrophone  channel.  The  decoded  ACODAC  amplifier  gains  are 
30,  to,  30,  20,  20,  and  20  dB  on  hydrophone  channels  1 through  6, 
respectively. 

E.  Shot  Data  Analysis 

1.  Introduction 


Shots  were  dropped  as  sources  at  ranges  of  10,  20,  and  100  miles 
from  the  two  AC0DACS.  For  each  range,  at  least  two  shots  each  were  used 
at  depths  of  300  and  QOO  ft.  The  actual  shot  pattern  as  determined  fro© 
the  ACODAC  tapes  is  shown  in  Table  VI. 


TABLE  VI  - SHOT  PATTERN 


Kang©  from 
Arrays,  miles 


10 

20 


50 


Number  of  Shots 
at  each  Depth 
300  ft  3-00  ft 
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For  digitizing,  the  analog  ACODAC  tape  was  passed  through  a 
10  to  300  Hz  bandpass  filter.  The  sampling  rate  of  600  Hz  was  obtained 
from  the  12th  harmonic  of  the  50  Hz  carrier  on  the  time  code  channel. 

This  phase  locked  sampling  technique  eliminated  wow  and  flutter  introduced 
by  the  recording  and  playback  processes.  For  each  shot  sequence  of  two 
to  four  shots  (at  a given  range),  the  digitizing  process  was  started  by 
using  the  first  shot  in  the  sequence  as  a trigger.  After  the  trigger, 

20  consecutive  blocks  of  8192  samples  (13.63  sec  at  the  600  Hz  sampling 
rate)  were  obtained  for  a total  continuous  digital  record  length  of 
273  sec.  An  example  of  one  such  sample  sequence  is  shown  in  Figs.  22 
and  23.  These  illustrate  the  signals  from  hydrophone  1 of  the  WH0I  ACODAC 
for  shots  at  3°0  ft  and  a 10  mile  range.  The  illustrations  are  from  hard 
copies  of  an  off-line  CUT  presentation  of  the  digitized  data.  Four  shots, 
including  the  trigger  shot,  are  clearly  shown.  Recording  system  gain 
changes  are  obvious  in  the  figures. 

2.  Multipath  Structure  of  Signals 

The  second,  third,  and  fourth  shots,  in  Figs.  22  and  23,  exhibit 
multipath  structure  with  at  least  three  large  stable  multipath  signals 
following  the  principal  arrival.  In  Fig.  2U,  an  expanded  presentation  of 
the  data  for  the  second  shot  is  shown.  Also  shown  is  the  output  of  the 
same  hydrophone  for  the  800  ft  shot  at  the  same  10  mile  range.  The 
multipath  structure  is  similar  for  the  two  depths  except  that  the  third 
and  fourth  arrivals  are  earlier  for  the  800  it  shot.  These  third  and 
fourth  arrivals  also  appear  to  consist  of  two  destructively  interfering 
signals  for  the  600  ft  shot. 

3.  Signal -to-Soise  Hat jo 


The  principal  arrival  for  the  second  shot  in  each  sequence  of 
shots  was  selected  for  S/ll  analysis.  The  S/H  was  estimated  according  to 
two  definitions  of  3/N  (Fig.  ?3).  The  two  S/'K  estimates  arc  computed 
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using  the  peak  of  the  signal  and  the  average  energy  of  the  signal. 
Based  on  the  previous  AKL  examination  of  the  shot  S/N  (Section  I),  an 
averaging  time  of  0,33  sec  was  used  in  estimating  the  average  noise 
intensity. 


In  Fig.  26,  the  amplitude  versus  time  plot  is  shown  for  the 
second  shot  of  each  sequence,  received  by  hydrophone  1 of  the  WHOI 
ACODAC.  Also  shown  in  Fig.  26  are  the  S/N's  for  each  of  the  illustrated 
shots.  The  temporal  limits  used  in  computing  average  signal  energy  are 
indicated  in  the  illustration.  This  same  set  of  information  is  given 
for  the  signals  from  the  remaining  hydrophones  in  Appendix  A. 

The  S/N  results  for  all  hydrophones  are  summarized  in  Table  VII. 
Some  clearly  discernible  trends  in  these  data  are  further  illustrated  in 
Figs.  27  through  3b. 

In  Fig.  27,  the  peak  S/N  is  shown  versus  hydrophone  number 
(WHOI  ACODAC)  for  the  10  mile  range  shots.  Subsequent  figures  show  these 
ratios  for  successively  greater  ranges:  20  miles  (Fig.  28),  50  miles 

(Fig.  29),  and  100  miles  (Fig.  30).  At  the  10  and  20  mile  ranges,  the 
near  surface  hydrophones  exhibit  the  best  values  of  S/N,  with  S/N 
decreasing  with  increasing  depth  of  the  hydrophone.  At  the  50  mile  range, 
a greatly  reduced  S/N  is  evident  in  the  data  from  hydrophone  3,  with 
values  from  the  deep  hydrophones  almost  equaling  the  two  shallowest  ones. 
At  the  100  mile  range,  hydrophones  2,  5,  and  b possess  the  higher  S/N 
values,  with  the  ratio  falling  off  both  above  and  below  these  depths. 
Clearly,  the  most  consistently  high  S/N  values  are  exhibited  by  hydro- 
phone 2,  though  hydrophone  1 shows  higher  S/N  for  some  ranges  and  shot 
depths . 


These  conclusions  are  further  illustrated  in  Fig.  31,  where  3/N 
values  for  hydrophones  1 and  2 are  plotted  versus  range.  The  greatly 
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Range/ Depth  of  Charge 
10  Biles,  300  ft 


TABLE  VII 


NAL-TO-NOISE  RATIOS:  SHOT  DATA 


Hydrophone 


W1 _ H2  W3  N4  KS  W6  M6 


27.8 

28.2 

28.7 

24.6 

20.6 

11.1 

21.3 

39.9 

37.9 

33.0 

32.1 

22. S 

28.7 

23.5 

KDU| 

mmm 

12.8 

if.  4 

32.3 

mSEk 

25.6 

27.9 

25.4 

2717 

23.9 

22.0 

mimm 

7.11 

39.2 

32.9 

38.1 

31.9 

23.9 

19.2 

26.8' 

22.4 

25.2 

23.7 

16.0 

14.8 

19.3 

33.3 

3S.3 

36.2 

26.8 

27.8 

28.8 

19.8 

16.8 

8.84 

15.7 

11.7 

■QQI 

S.IS 

29.2 

29.2 

19.6 

25.3 

24.3 

18.9 

iHM 

21.1 

£.16 

19.1 

21.6 

15.5 

36.2 

34.8 

14.5 

31.9 

35.2 

27.6 

13.9 
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FIGURE  29 
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FIGURE  31 
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decreased  S/N  for  hydrophones  3 and  4 is  evident  in  Fig.  32,  which  shows 
S/N  versus  range  for  these  two  hydrophones.  The  same  information  is  shown 
for  hydrophones  5 and  6 in  Fig.  33. 

In  Fig.  34,  S/S  values  are  compared  for  hydrophone  6 of  both 
the  WHOI  and  UM  ACODACs.  Although  comparable  S/N  was  exhibited  by  the 
two  systems  at  the  two  shorter  ranges,  at  the  long  ranges  the  S/N  was 
consistently  higher  at  the  WHOI  ACODAC,  especially  for  the  deeper  shots. 

4.  Spectrum 

The  intensity  spectrum  was  estimated  for  the  same  shots  used 
in  the  S/N  analysis  (the  second  shot  in  each  sequence).  In  Figs.  35 
and  36,  the  spectra  (2.34  Hz  resolution)  are  shown  for  300  ft  and  800  ft 
shots  received  by  hydrophone  1 at  the  10  mile  range.  The  same  information 
is  presented  for  other  ranges  and  hydrophones  in  Appendix  B.  Though 
there  is  considerable  variation  in  the  fine  structure  of  these  spectra, 
they  all  agree  in  general  form  by  having  a low  frequency  maximum  and  a 
general  decay  with  increasing  frequency  to  the  bandpass  filter  limit  of 
300  Hz. 
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Received  Signal  Amplitude  Versus 
Time  and  Signal-to-Noise  Ratio 
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APPENDIX  B 
Shot  Spectra 


Label  Example: 
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MEMORANDUM  FOR  DISTRIBUTION  LIST 

Subj : DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

Ref:  (a)  SECNAVINST  5510.36 

Enel:  (1)  List  of  DECLASSIFIED  LRAPP  Documents 

1 . In  accordance  with  reference  (a),  a declassification  review  has  been  conducted  on  a 
number  of  classified  LRAPP  documents. 

2.  The  LRAPP  documents  listed  in  enclosure  (1)  have  been  downgraded  to 
UNCLASSIFIED  and  have  been  approved  for  public  release.  These  documents  should 
be  remarked  as  follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  the  Chief  of  Naval 
Operations  (N772)  letter  N772A/6U875630,  20  January  2006. 

DISTRIBUTION  STATEMENT  A:  Approved  for  Public  Release;  Distribution  is 
unlimited. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 
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